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Abstract 

We investigate the splitting and mixing of p and uj mesons in nuclear matter. The 
calculations were performed on the basis of QCD sum rules and include all operators 
up to mass dimension-6 twist-4 and up to first order in the coupling constants. Special 
. attention is devoted to the impact of the scalar 4-quark condensates on both effects. In 

| nuclear matter the Landau damping governs the p — uj mass splitting while the scalar 

■ 4-quark condensates govern the strenght of individual mass shifts. A strong in-medium 

Q . mass splitting causes the disappearance of the p — uj mixing. 

1 Introduction 

m ; 

The investigation of in- medium modifications of hadrons is currently a topic of wide interest. 
. This is because the issue is related to chiral symmetry restoration as well as to a change 
of vacuum properties, and the phenomenon "mass of particles". Among the promising 
<3 candidates for a search for changed hadron properties in an ambient strongly interacting 

medium are vector mesons. Due to their decay mode V — > 7* — > e + e~ and the negligible 
interaction of the escaping e + e~ one can expect to probe directly the parent vector meson 
V. Indeed, strong evidences for changes of the p meson are found in relativistic heavy-ion 
collisions, where a meson-rich hot medium is transiently created (cf. 0). As the vector 
meson properties are coupled to various condensates [2 El El, which change as a function of 
both the baryon density and the temperature, complementary investigations of their behavior 
via the e + e~ decay channel in compressed nuclear matter is also required. Experimentally 
^ this will be done in a systematic way with the detector system HADES jS]. The situation is 

quite challenging since various predictions differ in details. 

In the invariant mass region up to 1 GeV there are various sources of e + e~ jUl El : Dalitz 
decays of many hadrons, bremsstrahlung, and the direct decays V — > e + e~ mentioned above. 
One important channel for di-electron production is the reaction tttt — > 7* — > e + e~. This 
channel has been evaluated with increasing sophistication over the last years (cf. pQ). The 
corresponding di-electron production rate R = dN ee /d 4 x in a medium characterized by the 
baryon density n and temperature T is given by (HI E] 

^iM^^^MiTKj^) (i-tg) . (1) 



dM ee K ee ' ' ' (2tt) 4 ee \ T J \ Ml 

Here K\ is a modified Bessel function, M ee stands for the invariant mass of the di-electron 
pair, and o is the total cross section of the process -k + -k~ — > 7* — > e + e~ 



4 aL, L Ami ,, . 



v{d\n) = -TT^f Jl - |F^,n)|^ , (2) 
6 q z y Q 

where F 7r (g 2 ) is the pion formfactor and q 2 = M e 2 e is the momentum squared of the decaying 
virtual photon 7*. 
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The p — uj mixing in vacuum has been discovered as a particular form of the pion formfactor 
l-F^g 2 )! a few decades ago jTU]- Since that time much work has been done aiming at a 
theoretical understanding of this mixing effect (for a review see [H]). Despite that there is 
still some debate concerning details of the p — to mixing in vacuum ^21 E3 El the experi- 
mental pion formfactor in vacuum can well be reproduced by means of several theoretical 
approaches. However, up to now the mixing effect has not been studied systematically in the 
medium. One may argue that the p — uj mixing is a tiny effect in evaluating the di-electron 
emission rate of warm nuclear matter. Moreover, inspired by [T^], it seems to be a generic 
effect of its own interest, which should be analyzed in a dense medium. This is one issue 
of the present paper. In addition, we are going to investigate the p — uj mass splitting and 
mixing effect simultaneously on the same footing, i.e. we use the same parameter set and 
the same approach for evaluating both effects. 

Refs. |TI3[ini[Hl showed that p and uj mesons experience, within the QCD sum rule approach, 
quite a different in- medium behavior. Even in zero- width approximation a large p — uj mass 
splitting was found when neglecting terms in the operator product expansion (OPE) which 
differ for p and uj mesons. The individual mass shifts depend on the yet poorly known 
density dependence of the four-quark condensate, assuming the same effective four-quark 
condensate for both the p and the uj mesons. A further goal of the present paper is to 
include all terms in the OPE up to mass dimension-6 and twist-4 (up to order a s ) and to 
study the importance of condensates which make p and uj differ. We extend our previous 
studies [TBI E| to consider here the yet unexplored effect of the density dependence of the 
four-quark condensate on the p — uj mixing in medium. 

Our paper is organized as follows. In section 2 we determine the mass shifts of p and uj 
mesons. We spell out the basic steps of QCD sum rules in low-density approximation and 
list all terms of the operator product expansion (OPE) up to mass dimension-6 and twist- 
4. We then present a numerical evaluation of the QCD sum rules and show that terms in 
the OPE which make p and uj differ are small at nuclear matter saturation density. The 
p — uj mass splitting is found to be determined by different Landau damping terms, while the 
individual mass shifts are governed by the density dependence of the four-quark condensate. 
The knowledge of the in-medium p, uj mass parameters is a prerequisite of a consistent 
treatment of the p — uj mixing studied in section 3. We define the phenomenology of the 
mixing and explain how this effect is related to observables. Afterwards we specify the QCD 
sum rule for the p — uj mixing and present details of the evaluation. The summary can be 
found in section 4. 

2 p — uj mass splitting 

The masses of p and uj mesons differ in vacuum by a small amount, Am = m w — m p = 11 
MeV. It was one success of the QCD sum rule method to explain this mass splitting in 
vacuum by differences in the OPE of p and uj current— current correlators [3]. Indeed, up 
to mass dimension-6 there is only one operator in vacuum, the so-called flavor mixing scalar 
operator, which differs in sign in the OPE of p and uj correlators. In the following we will 
investigate the behavior of this splitting at finite density, where nonscalar condensates play 
also a role. 
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2.1 QCD sum rule 



Within QCD sum rules the in-medium vector mesons V = p, oo are considered as resonances 
in the current— current correlation function 

ng?(g,n) = ijd'x e^ (0| T tf(0) |0) , (3) 

where g M = (go, q) is the meson four momentum, T denotes the time ordered product of 
the meson current operators J^(x), and \ stands for a state of the nuclear medium. In 
the following, we focus on the ground state of baryonic matter approximated by a Fermi 
gas with nucleon density n (in Ref. it was shown that temperature effects for T < 100 
MeV are of subleading order and may be neglected for our purposes). We first study isospin 
symmetric nuclear matter and extend later on our approach to asymmetric nuclear matter. 
In terms of quark field operators, the vector meson currents are given by 

J ^ = ^(ui^uTd'y^d) , (4) 

where the upper sign stands for the p meson while the lower sign stands for the u meson. 
We will keep this notation troughout the paper. Note that the interpolating currents eq. (@J) 
are based on the same field operators u, d. Therefore, evaluating the r.h.s. of eq. (j3J) will 
deliver the same condensates, however a few of them with different signs. To highlight this 
point we spell out all terms arising from eqs. (jSJ EJ in the following. 

We consider the nucleon and vector meson at rest, i.e. q^ = (go, 9 = 0) and = (Mat, k = 
0), which implies the vector meson to be off shell while the nucleon is on shell. Then the 
correlator (J3J) can be reduced to |lT^(g 2 ,n) = IT y )(g 2 ,n). In each of the vector meson 

V=p,U) 

channels the correlator U^ v \q 2 ,n) satisfies the twice subtracted dispersion relation, which 
can be written with Q 2 = —q 2 = —E 2 as 

U^(Q 2 ,n) _ n^)(0,n) ^w^l 7 . M^(^) (5) 

with Il( y )(0,n) = Il( y )(g 2 = 0,n) and nW(0) = dn( V ^f ,w) [ g 2 =0 as subtraction constants. 

We use IlW(0,n) = n/(AM N ) and 11^(0, n) = 9n/(AM N ) [TJ EE], respectively, which are 
the Thomson limit of the VN scattering process and correspond to Landau damping terms 

M- 

As usual in QCD sum rules |3J Hj, for large values of Q 2 one can evaluate the nonlocal 
operator of eq. (J3J) by OPE. We truncate the OPE beyond mass dimension-6 and twist-4 
and include all terms up to the first order in a s in the SU(2) flavor sector: 

T\(V)(Q2 n _ yAV) rj{V) U (V) n (V) (6) 

= 'M 1 + T c 'l)^{i)-M< + ^) (7) 
+\ (i + f^l) & (n|("K,®* + mdao|n> + ^^lf G2 l fi ) ( 8 ) 

-^a»i( n l (^7sA a n ^ 75 A a u + d W X a d d^ l5 X a d) \Q) (9) 



3 



±7ia s — (Q\ (u W X a u dYl 5 \ a d) \n) (10) 

-1^-^(01 (u^X a u u^\ a u + d 1 p\ a d d^X a d) \n) (11) 

-^a s -^(n\ (u lfl \ a u d^X a d) , (12) 
(m^(Q|m u «(7^G^|fi)+ml (^K^ G^d|fi)) , (13) 

n£U 2 = \^ f ^ q ^(Q\Sf(G-G au )\Q) (14) 

~ (I " ? Cf is) ' ^9"< n IST (wy^« + 5 7„M |n> , (15) 

iffl,^ = -^n/^gW^|ST (GSDrDxG^W) (16) 
3 + 1F° F W 1 qs^V^IST (u^pD u D\D a u + d lfl D u D x D a d) |0) , (17) 

= ±Ii^/V(%2sT (« 7M75 A a U d7,75A a d)|fi> (18) 

-g^e<7V W St (u 7 ^75A a M M 7 ,7 5 A a M + ^7 5 A a rf ^7 5 A a rf) |fi> (19) 

"^flVW St (li'y fl X a u (u lu \ a u + d lv \ a dj) |n> (20) 

-^gWk 2 st (S^d (u 7 ,A a M + d 7i ,A a d)) |0> (21) 

-^^gV W&ST (u [Dp, G ua ] + 7 a 75« + d [D,,, G ua \ + 7 Q 75 d) \ty (22) 

---^gV(«|St (m u + m d dL^d) , (23) 



where rif = 3 is the number of active flavors at a scale of 1 GeV, and Cp = (n 2 — l)/(2n c ) = 
4/3 with n c = 3 as number of colors; Op V = ~ [ 7ai , 7w ]_. The strong couplings are related by 
« s =£ 2 /( 47r )- 

The SU(3) color matrices are normalized as Tr (A a A 6 ) = 25 ab , the covariant derivative is de- 
fined as Dp = dp + igA a ^X a /2 and G 2 = Gp U G a ^ where G a is the gluon field strength tensor 
(G^ = G a ^\ a /2). The dual gluon field strength tensor is defined by Gp V = €p Vpa G a pa \ a j2. 
The OPE for scalar operators up to mass dimension-6 can be found in For the twist- 2 
condensates we have included all singlet operators with even parity up to order a s (nucleon 
matrix elements of operators with odd parity vanish). Their Wilson coefficients can be 
deduced from [213 • 1 

The Wilson coefficients of the twist-4 operators in lines (|18|) ... (|22|) are given in jS], and 
for the twist-4 operator in line (|23|) it can be deduced from |21j . where it has been found 
that the term (|23|) has some relevance for twist-4 effects of nucleon structure functions. The 
Wilson coefficient of an additional dimension-6 twist-4 operator, STg [Dp, G va ]_ 7 °g (for an 

1 Despite the fact that twist-2 non-singlet operators occur in the OPE for electromagnetic currents |2U|. 
they are absent in the OPE of eq. ©. Twist-2 non— singlet condensates contribute, however, to p — uj mixing. 
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estimate of this condensate see vanishes fll\ |2H]. 

We emphasize that the only difference between p and uj mesons in the truncated OPE consists 
in the terms in lines fjlOj) and (fTHj) . As mentioned above the term in line (I10J) is responsible 
for the p — uj mass splitting in vacuum; the term in line (|18|) vanishes in vacuum. It is 
now our goal to analyze the in-medium difference of p and u mesons stemming from the 
OPE side. Most terms in lines (J7J) ... (j22J) may be evaluated using standard techniques |2"4] . 
However, what remains to be considered are the flavor-mixing condensates in the lines (jlOj) . 
(H2J), (fTSj) . the mixed quark-gluon condensate in line (fH?j) . the pure gluonic condensates in 
the lines ()14)) . (|T6|) and the twist-4 condensate in line (J23|) . The QCD corrections to order 
a s of the twist-2 condensates in lines ()15|) and (fTTj) have not been taken into account in 
previous analyses. 

The chiral condensate and scalar gluon condensate have been discussed in some detail in [23] . 
Details for the scalar four-quark condensates in lines ••• dl) are given in Appendixes 
A and B, where also further notations are explained. The twist-2 quark condensates (lines 
(|15|) and (|17j0 and the gluonic twist-2 condensates (lines (|14|) and ()16Jl ) are explicitly given 
in Appendix C. The twist-4 condensates (lines (|18|) ... (J23)0 are listed in Appendix D. 
Performing a Borel transformation |3] of the dispersion relation eq. (jSj) with appropriate 
mass parameter M 2 and taking into account the OPE (JUJ) one gets the QCD sum rule 

noo („, „, - 1 /* e -./«* = Co M2 + 1 . (24) 

1=11 ' 

For nuclear matter we utilize the one-particle dilute gas approximation in order to evaluate 
all relevant condensates in the nuclear medium, i.e. 

(Q\0\Q) = (O) + ^-(N(k)\O\N(k)), (25) 
where the nucleon states are normalized by (N(k)\N(k )) = (27r) 3 2E^ 5(k — k') with 



Ek = yk + Mjy . The scalar dimension-4 and dimension-5 condensates are given by 

m u (Q\uu\Q) = m u (uu) + ^cr^ n , (26) 

(Q\^G 2 \Q) = (^G\-\m%u, (27) 

(fil^tZa-^G^Ifi) = A 2 (m) + iA 2 ^n, (28) 

2 m u 

where we have introduced the sigma term = m u (N(k)\uu\N(k)) /M^ and A 2 ~ 1 GeV 2 . 
The chiral d quark condensate follows by replacing the u quark by a d quark. The nucleon 
sigma term [2B] is 2a n = o~n + a N- 

Inserting the explicit expressions for all condensates (cf. Appendixes A ... D) one gets for 
the coefficients ci 5 2,3 in eq. (|2^|l 

- - M i+ v c 'z)> (29) 

ci = -^(rnl + m 2 d ), (30) 
c 2 = + Cf ~i) ( m u(uu) Q + m d (dd) + a N n 
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+ 



24 
1 

4 
112 



<^G 2 ) - °-M» N n 



Ti 



5 a s 
487 



_ —^±C F ) Af +d) M N n-^ ti/ — ^ Mjv n 



16 



7T 



(31) 



c 3 



7r a s k (gg)o 



1 + 



±36)Q§a s ._ 2 

J TV 



81 

5 n 
12 4 - 



k m q (qq) 
a N 1 



7? 



1 + 



67 



s C F —) At 
192/ 4 



M^77 + 



205 a. 



77/ A\M%n 



864 vr 

+1m w 77 (Jif 2 + ^ - (1 ± I)**, + 
where A^" -1- ^ = A% + with n = 2, 4, 2(qq) = (uu} + (<it/%, and 2m 5 = m u + m rf . 



(32) 



2.2 Evaluation 

We define a ratio of weighted moments 



/ ds ImlT y )(s,n) e" 



s/M 2 



m 2 v (n } M , s v ) = ^ 

/ ds lmU( v )(s,n) s~ l e- s / M ' 2 
o 



for which the desired sum rule follows by taking the ratio of eq. 
respect to 1/M 2 as 



(33) 



to its derivative with 



m v (n, M 2 , sy) 



c M 2 [l-(l + ^) e - 



S V /M 2 ] C2 C3_ 



M 2 



ili 4 



C 



2i\i (i 



n( y )(o,n) ' 

M 2 



(34) 



where we have identified the highlying (continuum) contributions as — InxQ^(s > Sy, n)/s — 
tccq (sy is the continuum threshold). The meaning of the parameter m v as normalized first 
moment of the spectral function Iinn^ becomes immediately clear in zero-width approxi- 



mation, — ImIF y )(s < sy,n) = nFyS(s 



77). 



* 2\ 
V ) 



from where my 



m 



v 



follows, eqs. (JSHl EU) 



are the corresponding generalizations for the case of finite width, in the spirit of a resonance 
+ continuum ansatz. The mass equation ()34j) is commonly used for describing m v in vacuum 
[3 HZl 123 123 E3 EH, at finite temperature jHj and at finite density gl E3 HB] and will be 
subject of our further considerations. 

The sum rule is reliable only in a Borel window M^ in < M 2 < M^ ax . If M 2 is too small the 
expansion eq. (|2*4^) breaks down. On the other side, if M 2 is too large the contribution of 
perturbative QCD terms completely dominate the sum rule. We adopt the following rules 



for determining the Borel window 



EH1 EU EHi] : The minimum Borel mass, M^ in , is 



determined such that the terms of order 0(1/M 6 ) on the OPE side contribute not more 
than 10%. The maximum Borel mass, M^ ax , is evaluated within zero-width approximation 
by requiring that the continuum part is not larger than the contribution of the resonance 
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part, i.e. 



1 




—m- 



(35) 



The parameter Fy can be evaluated by means of the QCD sum rule eq. ([24)1 . The obtained 
results for vacuum, F p = 0.0110 GeV 4 ,F w = 0.0117 GeV 4 , are in good agreement with the 
relations F p = m 4 / ' g 2 ^ = 0.0130 GeV 4 and F u = 9 m 4 /£ 7 = 0.0138 GeV 4 , respectively, 
which follow from the Vector Meson Dominance (VMD) [3*1 ITT| "3*""]. 

The threshold sy is determined by maximum flatness of my(n, M 2 , sy) as a function of M 2 . 
These requirements give a coupled system of equations for the five unknowns M^ in , M^ ax , 
Fy, sy, my. The final parameters Fy and my are averaged to get Borel mass independent 
quantities. For any parameter P this average is defined by 



For considering the mass parameter splitting effect there is no need to distinguish between 
isospin symmetric and isospin asymmetric nuclear matter since all operators in the OPE 
eq. © are isospin symmetric operators. Accordingly, in the following we study isospin 
symmetric nuclear matter. 

Twist-4 condensates have been estimated in where data of lepton-nucleon forward scat- 
tering amplitude has been used to fix the parameters K^, K 2 , and K\ d in eq. (f3*2*|) . The 
corresponding system of equations is under-determined and therefore various sets for these 
parameters can be obtained. We have investigated all six sets from for these parameters 
and find only very small changes of the results. In Fig. we show the results obtained 
with the parameter set given in Appendix D. Since in "TUJ [""""J "2H1 EZj a strong effect of the 
density dependence of the four-quark condensate was found we show here results for various 
possibilities, parameterized by k,n introduced in Appendix A, eq. ()87|) . The mass parameter 
of the p meson decreases with increasing density for all kn, while the u meson mass param- 
eter decreases only for sufficiently large Kjy. Other QCD sum rule analyses |U EH CH] have 
obtained also a decreasing p mass parameter. An increase of the u meson mass parameter 
has been found in [T51"""""], where the correct Landau damping term was implemented. 
The flavor mixing scalar operators (i.e. MJfy-, see Appendix B), while responsible for the 
mass splitting in vacuum, play only a minor rule in matter. That means, discarding the 
terms ~ Qq in eq. ()32|) yields curves which are nearly identical with those represented in 
Fig. ["3 The poorly known scalar four-quark condensate governs the strength of individual 
mass shifts, while the strong mass splitting in matter originates mainly from the Landau 
damping terms (0,n), which differ by a factor 9 for p and u |15| . The outcome of our 
study is that terms in the OPE, which cause a difference of p and u mesons, are small in 
matter since the mass splitting is mainly determined by the Landau damping terms. 




(36) 



In the following we will skip the average sign. 



2.3 Results 
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Figure 1: Mass parameter my of to meson (upper curves) and p meson (lower curves) as a 
function of the density for various values of the parameter n N ((a): = 1) (b): k-n = 2, 
(c): kn = 3, (d): = 4). The solid curves are for the full set of terms in eqs. (J2H1 - 
while for the dotted curves the twist-4 condensates are discarded, i.e. i^dud = 0- 





3 p — uj mixing 

First, we briefly describe the mixing scenario considered in the following. We follow the 
arguments given in jllj . The mixing can be accomplished by 

(37) 

where the subscript I denotes isospin-pure states, and e is the mixing parameter. The mixing 
formula (|37)l is quite general. Extending the mixed propagator approach described in ^T] to 
the case of finite density one can obtain the following relation between the complex mixing 
parameter e and the nondiagonal selfenergy 5 puJ (q 2 ,n) via (cf. [TT] for vacuum, cf. [TH] for 
matter) 

c( n ) = S pu (q 2 ,n) 

m t{ n ) ~ m2 P { n ) + « ImE tJ (g 2 , n) — i ImS p (g 2 , n) 
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The nondiagonal selfenergy S pLU (q 2 ,n), and therefore also the mixing parameter e, is directly 
related to the pion form factor, given by 



F w (q ,n 



2 \ 1 Q 9p1TTT 



g n q 2 -m 2 p (n) -2imE p (g 2 ,n) 

2 -i 

1 § u(q 2 ,n) ^ . (39) 

guy Q 2 — m t{ n ) — i hmE^g 2 , n) pui ' q 2 — m 2 — i Im£ p (g 2 , n) 

Since the main contribution of the second line stems from the region q 2 ~ m 2 , rr? w one usually 
approximates the nondiagonal selfenergy (^(g 2 , n) in the pion formfactor by its on-shell value 
at q 2 = m 2 = 0.5 (m 2 + m 2 ). 

The nondiagonal selfenergy consists of an electromagnetic part and a hadronic part, 
8 P u(jn 2 ,n) = <5™(m 2 ,n) + 5^{m 2 ,n). Both contributions can consistently be isolated in 
theoretical as well as experimental analyses. The electromagnetic part comes from the 
process p — > 7* — > u and can be evaluated analytically jHU]- In the following we are going to 
investigate the density dependence of 8 puJ {rn 2 ^n). 



3.1 QCD sum rule 

The basic object for the p — us mixing in matter is the mixed correlator 

n^(g,n) = ijd*x e^(T J>(x) J?(0)>„ , (40) 

with the isotriplet and isosinglet currents from eq. (|4*j). It is straightforward to recognize 
that 

n^(g,n) = i<T ntf{q,n) = IF(q,n) - U d (q,n) (41) 

with 

n q (g,n) = J Ae^(T5(x) 7M g(x) g(0)7 M g(0))„ . (42) 
This scalar function satisfies the twice subtracted dispersion relation 

™ = ™-n-(o, n) + Q *I/ ds ^. («, 

The subtraction constant n pa; (0, n) vanishes in vacuum JT] as well as in case of symmetric 
nuclear matter [T5] . For asymmetric nuclear matter, IP w (0,n) = —3 a np n/(4M^) with 
«np = ( n n — n p )/n [15J, where n n and ra p are the neutron and proton densities, respectively, 
and n = n n + n p . 

The other subtraction constant U puJ '(0,n) = dnP ^§ '"^ | g 2 =0 does not contribute to the sum 
rule after a Borel transformation. It is convenient |2E] to subtract the pure electromagnetic 
contribution p — > 7* — > u from hadronic and OPE sides of the dispersion relation ()43j) . In 
doing so we arrive at a new function, denoted by Il pw , which satisfies the same dispersion 
relation eq. (J4*3l . 

For large values of Q 2 one evaluates the l.h.s. of eq. (|4*3*|) by the OPE. Due to large can- 
cellations of the pure QCD terms according to eq. (jUJ) one has now to include also the 
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electromagnetic contributions to the OPE in contrast to eq. (JSj), where the electromagnetic 
terms are neglegible compared to the QCD terms. Accordingly, up to mass dimension-6 
twist-4, and up to first order in a s and a em the OPE is given by (for vacuum cf. [12 EH], for 



matter cf. [T5] ) 

11 W ) — iJ -scalar + ii c Z=4,r=2 + ii d=6,r=2 + ii d=6,T=4 + • • •> 

Scalar = «™ ^ " ^ " ™d) (45 

+\{l + ^C F -^^{n\{m u uu-m d dd)\n) (46 

1 a 1 - 

+ ^ — T^^K 4 m ^ u ~ m d dd)\Q) (47 

~\^s^i ({Sl\u llxl 5\ a u un 5 X a u\Q) - (Q\d^ 5 X a d dYl 5 X a d\Q)) (48 

-\™'7vi {^\ui^ a u u^\ a u\Q) - (n\dj„X a d dj p \ a d\n)) (49 

2 1 

--TTttemgJ (4(fi|w 7At75 w u^^u\Vt) - (fi|d 7/ , 75 d cf 7 M 75 d|0)J (50 

-^vra cm ^ (A{Q.\u ltl u u 7 p u\n) - (fi|3 7/1 d S 7 ^|fi)) (51 

+^TJ^6 H ^l m « ^V^N^) - ™2 (^K da, u G^d\Q)) (52 

+ey^^ Qm^fiK ua^ u F^u\n) + ^m 2 d (Q\m d da^F pu d\tt)^) , (53 

nLV=2 = -Q- v^i)^ 9 ^ 01 ^ M^-W)|fi> (54 

+^ ^ ^j^^lSf (4U7m^ - d^D u d) |fl> , (55 

K= & , T =2 = (3 + V Cf 3o) ^8^V(fi|St (u llM D„D x D a u - d llx D v D x D a d) \Q) (56 
+ — 7^7^^V(^|ST (4 u ltl D v D x D a u - d lfl D u D x D a d) \Q) , (57 

nLV=4 = -^gV(fi|^st(u 7At A a ^7,A a M - d 1 ,\ a dd lv X a d)\n) (58 

-^(Z'V W St (u 7M75 A a M U 7 , 75 A a M - d W X a d d lul ,\ a d) \Q) (59 

~^<f<fWg t St ^[^,G ra ] + 7 a 75M -3[^,G, a ] + 7 a 75C /) |0) (60 

-g^gV(n|ST (m u uI^LU - m rf di^ZW) (61 

-T^^gV^be ST (4 « 7M ™z 7i , W - d lt ,dd lv d) \Q) (62 
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2 1 
27 Q6' 



q p q u (£l\g e ST (4 wy^u wy v ^u - d^^d d~/ u ~/ 5 d) \Q) 



D cm F 



(63) 

n 5 d) \Q) . (64) 



V = e 2 /(47r) is the electromagnetic fine structure constant, F ua stands for the electromag- 
netic field strength tensor, and the dual electromagnetic field strength tensor is defined by 
Ffj, v = t^vpoF 90 ' ■ The covariant derivative of QED is defined as Z)^ m = <9 M + ieA^. The QED 
contributions may be deduced from the QCD terms by the replacements A a /2 — > 1 (which 
implies Cp — > 1) and g s — > e q (e q is the electric charge of quark q), respectively. Not all of 
the condensates given above have been taken into account in previous evaluations: the terms 
in the lines (J33E2E3I), the QCD corrections in the lines (j4*o1 the QED corrections 
given in the lines (|55jl and (|57|) . and all twist-4 contributions in lines ()58l - l64j) have not been 
considered yet. 

The isospin breaking of the scalar u and d quark condensates is usually parameterized by 



7 + 1 



(0|dd|0) _ (N\dd\N) _ (Q\dd\Q) 
(0\m\0) ~ (N\uu\N) ~ (n\uu\Q) 



(65) 



where we have generalized the corresponding relation for vacuum |15| IT2"1 13*%] to the case of 
nuclear matter. 

The four-quark condensates are given in Appendix A. The twist-2 quark condensates are 
listed in Appendix C and the corresponding parameters can be found in Appendix E. The 
twist-4 condensates, listed here for the sake of completeness, are neglected in our analysis 
since they are strongly suppressed in the chosen Borel window. 
Performing a Borel transformation of eq. ()43j) leads to 



oo 

U puj (0,n) - - fds 

7T J 



lmU puJ (s,n) 



-s/M 2 



do M 2 + £. 



di 



i=i 



(66) 



The coefficients c?i,2,3 in linear density approximation and neglecting all twist-4 condensates 
are given by 



d 
di 

d-2 



;OL P 



64tt 3 em ' 
3 , o 



m 



8tt 2 



a_s_ 

71 



d) i 
1 

4 



(m u - m d ) + 



1 C^en 

72~ 



(Am u - m d ) 



(??)( 



n 288 



M N n [AV + A 



id,p 



+ df 



14 



-txko Vtqa a - a e 



(qq) 



1 , Kn 
1 H n 



a N 



k m q (qq} 0/ 



67 



7i 1152 



M 3 N n{A u 4 ' p + A 4 



d.p 



+ d 



AS 



2m„ 



-n 



(67) 
(68) 

(69) 

(70) 



where further terms proportional to m q 7, 7 2 and 7 a em have been neglected. The terms 
df s (I = 2,3) are proportional to a np and account for isospin asymmetric matter. Their 
impact on mixing will be considered separately in subsection 3.3.3. 
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Finally, we specify the hadronic side of the QCD sum rule eq. ([Sfljl (cf. [31121112] for vacuum, 
[To] for matter), where the <fi meson has been implemented in accordance with ^3] 



1 ImlP^n) 1 r 

71 s 4 

1 

+ 4 



fp S(s - m 2 p ) - f u 8(s - ml) + ft 6(s - m 2 ] 



a,, 



fa 5{s - mj) - fa 5(s - mi,)\ + ^ Q(s - s v ) . (71) 

The necessity for including the higher resonances p' and u' is discussed below. 
We mention that the term is allowed since the meson is not a pure ss state but mixed 
with the lo meson. Even more, it has been found in ^3] that the <fi meson gives a significant 
contribution in vacuum due to large cancellations between f p and fa. Accordingly, we drop 
the assumption of ideal mixing and take into account such a term. 2 

The five parameters J ' p , j ' p > , fa, fat and fa have to be evaluated selfconsistently within the 
QCD sum rule approach. What we still need is a connection between these new parameters 
and the parameter S pLU (m 2 ) which enters physical observables like the pion formfactor in 
eq. ([59]). Such a relationship can be obtained by means of VMD [TTj 

Tfl 2 TO 2 

J>(s) = -^<(x) , J${x) = 3 -^(s) , (72) 

where (fY(x) is the field operator of the respective vector meson V — p,ui. If one inserts 
these relations into the correlator eq. ([4Dj) one gets an expression which relates lifjf with the 
mixed propagator (keeping in mind the zero- width approximation at all stages). Another 
expression for H?" can be obtained by inserting eq. ([7T|) into the dispersion relation eq. ()43|). 
Equating both expressions leads to the searched relation 

O^ 2 ) = -(//» + /«) ^ </p7&>7 > ( 73 ) 

which is valid to order 0(Am 4 /m 4 ) being a fairly well approximation even when taking into 
account the strong mass splitting between p and ui mesons. We mention that for evaluating 
the momentum dependence of ^^(q 2 ) the applicability of VMD has been debated in [T3J due 
to the impact of the meson. On the other side, the reliablility of VMD for a momentum 
dependence of ^^{q 2 ) has been confirmed in |14j . where the finite width of vector mesons 
is taken into account. Anyhow, in zero- width approximation VMD, which leads to eq. ()73[). 
is applicable as long as one restricts oneself to the on-shell value of this quantity, i.e. to 
(m 



3.2 Evaluation 

Since the most relevant parameter S pui enters the approach via the combination £ ~ f p + fa 
it is convinient to rewrite the sum rule eq. (jMj) as 

1 m 2 _/m 2 __ \ 2/M2 1 rnT (™r_ g ,\ '/m* , IJ_ f 
eM 2 \M 2 P + ^M 2 \m 2 P + AM 2 ' 




OO J 

em - S y/M 2 - * ■* "' 



2 In finite width QCD sum rule, which is necessary when considering the momentum dependence of mixing 
fipuil 2 ), the 4> contribution is negligible JTJ]. We mention also that there is no need to use unphysical values 
for m p ' and m u i, as has been pointed out in 
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where we have introduced [T3] 



m' 



C m 4 l ' 2 ) ' C '" 7"7 ; l 2 



a _ /hj /p m2 at _ fp 1 m ' , 



with 2m 2 = (m 2 + m 2 ), 2m' = (m 2 , + m 2 ,), Am 2 = m 2 — m 2 and Am' 2 = m 2 , — m 2 /, 
respectively. We stress that eq. (|73J) is valid to order C(Am 4 /M 4 ). Despite the observed 
strong mass splitting found in the previous section, eq. (|71jl is a good approximation: The 
terms of order 0(Am 4 /M 4 ) would give less than 10 percent correction to terms of order 
0(Am 2 /M 2 ), even at such a small Borel mass like M ~ 1 GeV. 

The residues in the hadronic model (fTTj) can be expressed by the new variables (J75j) to give 

Finally we give an expression for the mixing parameter e in zero-width approximation (i.e., 
Im S PjW = 0) which can be deduced from (|T3|) and (}38j) 

2 

6 = "A^IT C • (78) 

We need five equations for the five unknowns (, ft, /3', fy. One could perform a Taylor 
expansion of eq. ()66|) ending up with an equation system for these five parameters. This is 
the frame work of Finite Energy Sum Rules (FESR). Instead, here we use a combined FESR 
and Borel analysis, following the approach described in jT21 El which we extend to finite 
density. Accordingly, the first equation comes from a local duality relation [TJ] which results 
into 

411^(0, n) - ft (m 2 - ft' ('ri 2 = 4d s v + 4d 1 - f<p (79) 

and agrees with the first equation of the FESR approach ^21 ITTTj . This equation makes 
clear why the higher resonances p' and u' have to be taken into account: Without these 
higher resonances one would get either ft w or Q ~ which would be in contradiction with 
experimental findings. The second equation is just the sum rule eq. Two equations are 
obtained by the first and second derivatives with respect to 1/M 2 of eq. (JBBJ), cf. [12] (due to 
the high Borel mass and the small contribution of the threshold term the second derivative 
sum rule is applicable, in contrast to the mass splitting, investigated in the previous section, 
where a second derivative sum rule becomes unstable (H3|). 

For evaluating we still need a fifth equation. In [T3j a third derivative of sum rule has 
been used which could cause instabilities due to the truncation of OPE (SHESj- To avoid 
such instabilities the individual contributions of p' and u' have been approximated by an 
effective strenght f p > w > at the averaged mass of m p i^i in [T3]. However, at finite density ft is 
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density dependent. Therefore, in order to improve this approximation we apply the second 
FESR for the parameter (3, i.e., 

(1 + (3) C m 4 + (1 + (3') C ^ 4 ~Uml = -2 d s 2 v + 4 d 2 . (80) 

The resulting system of equations has to be solved selfconsistently giving the five unknowns 
as function of the Borel mass, ((M 2 ), C'(M 2 ), (3(M 2 ), (3'(M 2 ) and U(M 2 ). 
In Fig. |21 we have plotted these parameters as a function of the Borel mass for different 
densities. Like in the Borel analysis for the p — uo mass splitting we have to find an 
appropriate Borel window M^ in , M^ ax . To determine the minimal Borel window one could 
use again the 10 % rule getting M min « 1 GeV, while in [T2| a 25 % rule has been used 
getting M min 1.3 GeV. But it turns out that in such a region around M min the sum rule is 
unstable for a wide range of parameters [T2| IT3j. Nevertheless, the curves in Fig. |2] evidence 
that a stable region for all five unknowns exists in the interval 4 < M < 8 GeV. This 
observation confirms a corresponding stability investigation in Therefore, in line with 
13 , we will use a static Borel window M m j n = 4 < M < M max = 8 GeV over which we have 
to average (using eq. J3HJ)) to get Borel mass independent quantities. The result found in 
f21E3E] that the threshold parameter so in vacuum turns out to play a subdominant rule 
is also valid in case of finite density. Accordingly, we may use a fixed value, sy = 2.0 GeV, 
for all densities. 




02468 10 02468 10 

M [GeV] M [GeV] 



Figure 2: Parameters ( (a), f3 (b), £' (c), (3' (d) and (e) as a function of the Borel 
mass. Dotted curves are for vacuum, solid curves stand for n = and dashed curves depict 
n = 2 no- All curves are evaluated for kn = kq = 3 (here the mass shift of vector mesons 
has not been taken into account). 
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3.3 Results 
3.3.1 Vacuum 

First of all, let us briefly discuss the pion formf actor in vacuum, given by (fHHj) with n = 0. 
Our sum rule analysis for vacuum results in ( = 1.055 x 10~ 3 in good agreement with 
Using (|75j) gives for the mixing parameter e = —0.21. The hadronic contribution of 
the nondiagonal on-shell selfenergy which enters the pion formfactor, is given, via eq. ()73|) . 
by S^(m 2 ) = — m 2 g p ^ g wl C/12 = —4289 MeV 2 which amounts, by taking into account 
the electromagnetic nondiagonal on-shell selfenergy S™(fn 2 ) = 610 MeV 2 |3*9*1 I12j. in total 
to 5 pU )(jn 2 ) = —3679 MeV 2 , in fair agreement with experiment [10] • Using this value we 
get the pion formfactor in vacuum as shown in Fig. El It reproduces very well recently 
obtained experimental data [IT]. In ^3] it was argued that the sum rules might give not 
a good agreement with data when taking the parameter set of [12]. Our analysis shows, 
however, that the sum rules are in agreement with experimental data when using appropriate 
parameters. 

50 
40 

CM 

— 30 
\£ 20 
10 


600 700 800 900 1000 
q [MeV] 

Figure 3: Comparision of the formfactor evaluated within the QCD sum rule method (solid 
curve) and the results of the CMD-2 experiment (symbols) [4*T] . 




3.3.2 Isospin symmetric nuclear matter 

After reproducing the pion formfactor in vacuum we now turn to the density dependence 
of the mixing effect. Due to the small effect of mixing compared to splitting and the large 
impact of the four-quark condensate and Landau damping terms on mass parameter splitting, 
it becomes obvious that mixing does not strongly influence the mass parameter splitting 
effect. But on the other side, the mass parameter splitting effect could strongly influence the 
mixing effect. To study the effect of the p — uo mass parameter splitting on the p — io mixing 
we have to implement in the five equations for the five unknowns (, (3, /3', the density 
dependent mass parameters, i.e.: m p (n), m^n) and m^n), respectively. For m p (n), m^n) 
we use the values obtained in the previous section, while for the density dependence of the 
(f) meson we will take the relation m^n) = (1 — a n/no)m t f ) (0) with a = 0.03, which turns 
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out to be almost independent of kn ^Hl- The results for the five parameters are shown in 
Fig. H 




Figure 4: Parameters ( (a), j3 (b), (' (c), (3' (d) and (e) at finite density. Dotted curves 
are for = 2, dashed curves are for = 3 and solid curves are for = 4. The density 
dependence of the mass parameters of p, uo and mesons (without the twist-4 condensates) 
has been taken into account consistently. 

From the density behavior of the parameter ( (see Fig. 0] (a)) one might conclude that the 
mixing effect remains in matter. But this is actually not the case. In view of eq. (J78j) we 
recognize that the mixing angle e is strongly suppressed by the factor 1/ Am 2 . Additionally, 
the mass shift of the p meson modifies significantly the pion formfactor. Using eqs. (139(1 and 
(JTJ) for the pion formfactor and di-electron production rate, respectively, we get the results 
shown in Fig. 

These figures show that the mixing effect in the pion formfactor as well as in the di-electron 
production rate is washed out due to the mass shifts of the vector mesons. But one has 
to keep in mind that global changes of vector mesons in matter like mass shift and width 
broadening turn out to be correlated in nuclear matter [321 EHJ EH] • Taking into account 
such broadening effects needs further investigations. 

We also show results without the mass shifts of p, u and mesons. The corresponding 
density dependence of the five parameters £, /3, £ , (3 , is shown in Fig. One observes 
noticeable changes for the parameters. The dashed curve (i.e. = n ) in Fig. El (a) recovers 
the density-independence of ( for isospin symmetric nuclear matter as anticipated in [To] . 
Otherwise, depending on the parameter which governs the density dependence of the 
four-quark condensate, ( may slightly increase (large k^) or decrease (smaller k^) with 
increasing density. The resulting pion formfactor and the di-electron production rate are 
plotted in Fig. At finite density one obtains a very small modification of the formfactor 
compared to the vacuum, while the modification of the rate is nearly invisible. 
In the Figs. |S] and [7| the meson peaks are assumed to be distributed with a schematic 
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Figure 5: Left pannel (a): Pion formfactor at saturation density n = uq. Mass shifts of 
vector mesons (without twist-4 condensates) are taken into account. The dotted curve is for 
vacuum, while the solid curves are for saturation density n = n$. The labels 1,2,3 denote 
k n = 1,2,3, respectively. Right pannel (b): Di-electron production rate from pion-pion 
annihilation at finite density and T = 100 MeV. The dotted curve is for a hot pion gas and 
baryonic vacuum n = 0, while the solid curves are for saturation density n = n . 
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Figure 6: Parameter ( (a), (3 (b), (' (c), (3' (d) and (e) at finite density. The density 
dependence of the mass parameters has not been taken into account. Same notation as in 

Fig. a 



width JmE p (E) = -^je(E 2 - 4ml) 3 / 2 @(E - 2m w ) and ImS^E) = -m^T^ Q(E - 3m,), 
respectively. In Fig. [o] the density dependence of taken into account, while in Fig.[7| 

no shifts of assumed. Obviously, the di-electron rates shown in Figs. |S] and [7| differ 

significantly. 
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Figure 7: Left pannel (a): Formfactor at finite density. Dotted line denotes vacuum, dashed 
line represents n = no and solid line means n = 2uq. Right pannel (b): Di-lepton production 
rate for pion-pion annihilation at finite density for T = 100 MeV. The plotted curves are for 
ktv = 3. No mass shifts 



There is the possibility, advocated in [32], that in- medium the original (vacuum) p peak is not 
shifted, but additional strengths develops below the p peak. A similar possibility has been 
reported in [7] for the u) meson. In such cases the p — uj mixing remains, similar to Fig. [3 but 
the weighted p strength is shifted down, as required by the sum rule considered in section 2. 
A proper handling of this situation deserves further investigations with explicit knowledge 
of the p and u in medium spectral functions. Experimentally, precision measurements with 
HADES 5] can deliver informations on the in-medium behavior of the p — uo mixing. 



3.3.3 Isospin asymmetric nuclear matter 

So far we have considered isospin symmetric nuclear matter. While it is not necessary to 
study isospin asymmetric matter for the mass splitting effect, finite values of a np have some 
relevance for the mixing effect Therefore, in this subsection we concentrate on isospin 
asymmetric nuclear matter. The needed proton and neutron condensates are given in the 
Appendix E. Accordingly, the coefficients df s in lines (I69|) and (|7U|) contain the following 
terms proportional to a np . 

4 S = ~\ (l + ^ c f 7^ m <i a ^P (P\uu - dd\p) 

~h ~V W~ \ anp 5 mq ^ U ~ 



1 5 



| |s _ C F ) [AY-Ai P ) M N a np n 
+ — §7M N na np (Ar-A d 2 ' p ) , (81) 

d 3 S = 37 71 a * a n P n -rjr (99)o ~ dd\p) 
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7 1 

+ n a ™P n TT~ Kn (Wo 0°I mm - dd\p) 

+7^na„ p (Ar-A d f) , (82) 

where terms of order O ((m^ — m u )a np ), O (7a em ) and O {pf a np ) are neglected. 

The dependence of the parameter (, eq. (J75J), which governs the mixing effect of the pion 

formfactor via (J73J), on the asymmetry parameter a np is seen in Fig. |H] (a). For strong 

asymmetry one obtains a remarkable increase of roughly linear with a np . We note that the 

dashed curve of Fig. |H] (a) is in good agreement with ^3] where an asymmetry dependence 

( = ((°) + a np n/(0.2n ) with C (0) = 1-1 x 10~ 3 and C (1) = 1-5 x 10~ 3 has been reported, 

while our findings correspond to = 1.05 x 10~ 3 and = 1.9 x 10~ 3 . 

Altogether, without accounting for the mass shifts, an amplification of the mixing effect in 

the pion formfactor is obtained (see dashed curve in Fig.|H](b)). In contrast, when accounting 

for the individual mass parameter shifts the mixing effect is washed out (see solid curve of 

Fig. El (b)). 




Figure 8: Left pannel (a): Parameter ( as a function of a np at saturation density Uq (solid 
line: individual mass shifts of vector mesons have been taken into account; dashes line: 
without mass shifts of vector mesons). Right pannel (b): Pion formfactor for = 3 and 
ot np = 0.2 (dotted line: vacuum; solid line: n = n with mass shifts of vector mesons; dashed 
line: n = n without mass shifts of vector mesons). 

Finally, it is expedient to summarize the differences between the analysis presented here and 
Ref. |15j . which are, so far, the only investigations where the QCD sum rule approach has 
been applied to the mixing effect at finite density. Besides the usage of a complete OPE up 
to mass-dimension-6 twist-2 for the mixing effect and a selfconsistent Borel analysis for all 
unknowns at finite density in our work, the improvements are the following: First, we have 
implemented the individual mass parameter shifts of the vector mesons in a consistent way 
and have studied their impact on the mixing effect. A second difference consists in taking 
into account the meson on the hadronic side, which is necessary due to large cancellations 
between p and u> meson contributions (this has been pointed out for vacuum in ^H])- Thirdly, 
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we have investigated the relevance of ( by considering the influence of the mixing on pion 
formfactor and di-electron production rate. 

4 Summary 

In summary, we have investigated the mass parameter splitting and the mixing of p and uj 
mesons in nuclear matter within the QCD sum rule approach, starting from a complete OPE 
of the current-current correlator up to mass dimension-6 twist-4 and up to the first order in 
the coupling constant. Special attention is devoted to the impact of the poorly known scalar 
4-quark condensates. We have found a strong p — uj mass parameter splitting. The scalar 
flavor mixing condensate has been evaluated at finite density using quite general assumptions. 
It turns out that this condensate, while responsible for the p — uj mass parameter splitting 
in vacuum, plays a subdominant role in matter. Instead, the individual mass parameter 
splitting of p and uj mesons is mainly governed by the Landau damping terms. The scalar 
4-quark condensates have a strong impact on the individual strengths of the mass parameter 
shifts, while the amount of the splitting is fairly insensitive to these condensates. 
We emphasize that the mass parameters are weighted moments of the spectral functions. 
A mass parameter shift in medium does not necessarily mean a simple shift of the peak 
position of a spectral function, rather additional strength may occur at lower or higher 
energies causing a shift of the weighted moment. The presently employed form of the QCD 
sum rule approach is not sensitive to such details. Only a detailed modelling of the hadronic 
in-medium spectral function with parametric dependences allows for more concise statements 

Another physical effect investigated concerns the p — uj mixing at finite density and the 
impact of the p — uj mass parameter splitting. Starting with the vacuum we find an excellent 
agreement with experimental data recently obtained. In medium, the nondiagonal selfenergy 
Sp^im, n), which drives the mixing effect, is only weakly amplified in isospin symmetric 
nuclear matter. The mixing parameter (, however, is remarkably enlarged for strongly 
isospin asymmetric nuclear matter, such as in uranium nuclei with a np = 0.2. Therefore, 
not taking into account the individual mass shifts of the p and uj meson would indeed result 
in an in-medium amplification of the mixing effect. In contrast, if one takes into account 
the strong mass parameter splitting of p and uj mesons as a pronounced splitting of the 
corresponding peaks then the mixing effect in the pion formfactor as well as in the di-electron 
production rate disappears in medium, both for isospin symmetric and isospin asymmetric 
nuclear matter. Upcoming measurements at HADES can deliver valuable information on 
these issues. 
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A Scalar flavor- unmixing 4-quark condensates 



In lines (JSJ) and (jTTjl one recognizes two different types of scalar flavor-unmixing 4-quark 
condensates (q = u, d) 



Mf = {n\q lfll ,\ a qqYl,\ a q\£l) 



(83) 



and 



M™ = (n\q^\ a qqY\ a q\ty . 



(84) 



Previous studies employed a factorization for the scalar flavor-unmixing 4-quark conden- 
sates |2E]- We go beyond such approximation, = yft(fi|gg|f2) 2 , pointing out that k 
is uncertain and might even have a density dependence. In the spirit of the linear density 
approximation eq. a Taylor expansion results in 



M q A « = j(qq)l4 



1 + 



(i) 

k n ' a N n 
m q (qq} 



(85) 



The first term, i.e. y^g) 2 , k o\ * s mere ly an expression for {q^^^,\ a q qj ,1 'ysX a q}o. The 
second term, proportional to , parameterizes the poorly known 4-quark condensate in 
the nucleon (N(k) 1^7^75 X a u u r y' 1 'y5X a u\N(k)). Similarly, for the other 4-quark condensate 
we obtain 



My 



16- x2 (2) 



(2) 

k n ' a N n 

4 2) m i (qq)o 



(86) 



Accumulating all flavor-unmixing four-quark condensates, with the right weight given from 
the OPE, one obtains in linear density approximation finally 



-MT - -Mf - -M v u - -Mt d 
2 A 2 A 9 v 9 



M 



7 16, > 2 



1 + 



km o n n 
m g (qq)o 



(87) 



Note that kjv = Ko is conform to the large- N c limit [13] • Since we are interested in medium 
effects, we adjust the value of to the vacuum masses, yielding k = 3 both for p and u, and 
study the impact of the unknown parameter k n . As stressed in [111 EI], only a comparison 
with experimental data can pin down k^. 
For treating the p — u mixing we also need 



Nf 



(88) 



and 



K q 



{n\q-f^q qYq\ty . 
With the same steps as above we arrive at 



(89) 



4A^ U + N A 



dd 



9 



+ - Ny d 



(qq} K o 



K N (J N n 

«o m q (qq)o 



(90) 
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B Scalar flavor- mixing 4-quark condensates 

Now we estimate the two scalar flavor-mixing condensates in lines (|10jl and f!12[) at finite 
density. (For the vacuum such an estimate is given in [2E]-) Let us first consider the 
condensate in line (fTUj) . To evaluate such a condensate we insert a complete set of QCD 
eigenstates after a Fierz transformation 

Mf = ({i\u W \ a ud^ l5 \ a d\n) 

= £<n|n?df|n)<n|2J^|n) (A 8 )« 

n 

x (fftfrP + i( 7/1 )^( 7 ^)^ + I( 7m75 )^( 7 m 75 )7/3 _ ( 75 )^( 75 )-r^ , (91) 
and approximate the sum by 

£ |n)H « |o><o| + |n*><n*| + E/ 7^ ^ 1^ Ap))<Ap) "I , (92) 



where is the ground state of matter, denotes low-lying excitations (e.g. particle-hole 
excitations), and \Q7c b ) means ground state plus pion with isospin index b (other states with 
mesons heavier than pions are suppressed by their larger masses). The matrix elements 
(Q\v^di\fl) and (n|uf df vanish due to quark flavor conservation yielding 

x (yapp + 1 (7m) «^ M 7/J + 1 ( 7m75 )^( 7 ^ 75 )^ - ( 75 ) a5 ( 75 )^) . (93) 

The soft pion theorem [HJ ^] allows to calculate the needed terms in linear density approx- 
imation as 

<n|B?df |n tt^o)) = ^^(75) Q5 (fi|gg|fi) , (94) 

(fi|Uf4|fi7r 2 (0)) = -ll^( 75 )^(fi|gg|fi) (95) 

and (Q\v%di\fl tt 3 (0)) = 0. Inserting these matrix elements into (f9"3|) results in 

i< = (n\qq\n) 2 Ql (96) 

with the cutoff Qq coming from the momentum integral. With eq. (J2*5*|) and (N(k) \qq\N(k)) = 
one gets the final result 



T ud _ 4 lg n \ 2 



1 0"V^ 



(97) 



to 9 <9?)o. 

Using the same technique for the matrix element in line ()12|) one finds 

= (n\u lf ,X a u d^X a d\n) = -Mf . (98) 

Qq is adjusted to the vacuum p — u mass splitting. Using Q = 150 MeV we get the right 
experimental vacuum masses, i.e. m p (0) = 771 MeV and 771^(0) = 782 MeV [13] for our 
chosen parameters, (qq) = (—0.245 GeV) 3 , (a s /ir G 2 ) = (0.33 GeV) 4 , a s = 0.38, m u = 4 
MeV, m d = 7 MeV, M% = 770 MeV, a N = 45 MeV. Furthermore, we take the known 
vacuum values of Mat, /tt and m^. 
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C Twist-2 condensates 

The quark twist-2 condensates appear in lines (fT5|) and (fTTjl . respectively, while the gluonic 
twist-2 condensates appear in lines (fUj) and (|Tfi|) . respectively. The operator ST creates 
a symmetric and traceless expression with respect to the Lorentz indices, i.e., for spin- 
2 ST {O a p) = ^ ■ (O a p + Op a ) — \ g a pOJ 1 and analogously for spin-4 condensates. These 
condensates vanish in vacuum and therefore, according to the low-density approximation 
we need only the nucleon matrix elements which can generally be written as 



(N{k)\$Tq^D„q\N(k)) = -iS^ A q 2 (Li 2 ) , (99) 

(N(k)\STGfG a „\N(k)) = S, v A^ll 2 ) (100) 
for spin-2 operators and 

(jV(fc)|STg 7 ^D A £U|iV(fc)> = iS^ A\(^) , (101) 

(N(k)\SfG^D u D x G pa \N(k)) = -S^ Xa A°(p?) (102) 
for spin-4 operators, respectively. The Lorentz structures are defined as 

Sfiu = k^k u —k g pv , (103) 



s. 



fll/Xcr 



k 4 

k^kvkxK + — (g^gxa + g^g™ + g^g v \) 
k 2 

{k fl k u g\ a + kpk\g va + k^k^gxu + k u kxg fJi<J + kyk^g^x + kxk a g pv ) 



.(104) 



The reduced matrix elements of quark twist-2 condensates are defined as 
i 

A 9 (li 2 ) = 2 J dx x 1 ' 1 [q N (x, li 2 ) + q N (x, li 2 )}, where g7v(x, /i 2 ) and q N (x, li 2 ) are the 

o 

quark and antiquark distribution function inside the nucleon. We take A^^lGeV 2 ) = 

1.02 and A^ +d \lGeV 2 ) = 0.12 respectively. The reduced matrix elements of gluon 

i 

twist-2 condensates are defined by Af (/i 2 ) = 2 / dx x 1 ~ 1 Gn(x, li 2 ), with Gn( gluon 



distribution function inside the nucleon at the scale fx 2 . We use ^(lGeV 2 ) = 0.83 and 
^(lGeV 2 ) = 0.04 [21, respectively. 



D Twist-4 condensates 



Twist-4 condensates appear in the lines (jl8imil2Tj) and (J22j) . All twist-4 operators vanish in 
vacuum and therefore, according to the low-density approximation (J25|) . one needs only the 
nucleon matrix elements. The nucleon matrix elements of symmetric and traceless twist-4 
operators can be decomposed as 



{N{k)\ig s ST hi 



n n 

^ [ii ^ va 



7 a l5 u)\N(k)) 



(105) 
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(iV(fc)|^St(ll7 M 75A a wll 7l , 75 A a M )|iV(fc)) = 2S, v (Kl + K 1 d -K 1 ud ) , (106) 
(iV(fc)|^St(ll7^75A a nrf7,75A a rf)|iV(fc)) = 2 5^^), (107) 

(iV(fc)|^ 2 Sf (u^X a u(u^X a u + d^\ a d))\N(k)) = 2S, V (K 2 U + K 2 d ) (108) 

with S^ u defined in eq. ()103|) . The other twist-4 condensates, where u and d are interchanged, 
are equal to the given ones due to the assumed flavor symmetry. As pointed out in [HE] the 

1 2 

coefficients K u ' d ud are related to the nucleon forward scattering amplitude of the electromag- 
netic current. We take the following parameter set: K\ = —0.112 GeV 2 , K 2 = 0.110 GeV 2 , 
Kl = -0.300 GeV 2 , K x ud = -0.084 GeV 2 as default. For the d quark we use K 1 /' 9 = f3K^ 
with p = 0.476 from [HE]. 

We remark that the parameters ' 2 d 9 ud should be taken at a hadronic scale of p = 1 GeV. 
Unfortunately, twist-4 condensates are poorly known and even available only at a scale of 
p = 2.25 GeV. To evolve these parameters down to p = 1 GeV would require the knowledge of 
anomalous dimensions which are not available. Here we use the above condensates, expressed 
by K^, K^, and K^ d , to demonstrate that accounting these condensates has little influence 
on the mass splitting and individual mass shifts of p and uj mesons. 
For the twist-4 operator in line f!23l) we use the estimate [2*T] : 

(N(k)\m q qD^D v q\N(k)) ~ -P^(N{k)\m q qq\N(k)} , (109) 

where P^ is the average momentum carried by the quark q inside the nucleon. Taking 
Pjl ~ [21\ (k^ is the momentum of nucleon) and making the operator symmetric and 

traceless we get 

(N(k)\Sf m q qD^D u q\N(k)) ~ S^m q (N(k)\qq\N(k)} . (110) 

E Parameters for p — uj mixing 

For the p — uj mixing we have to distinguish between proton and neutron matrix elements. 
In particular we need the following twist-2 and twist-4 condensates 

(p(fc)|ST q llx D v q\p{k)) = -iS, u A q f , (111) 

(p(fc)|ST u ltl D u D x D a u\ P (k)) = iS, uXa Ar (112) 

with the proton state \p(k)) and analog expressions for the neutron. The reduced matrix 

i 

elements Af p (p 2 ) are defined as Af p (p 2 ) = 2 / dx x l ~ l [q p {x, p 2 ) + (— l) 1 q p (x,p 2 )}, where 

o 

q p (x,p 2 ) and q p (x,p 2 ) are the quark and antiquark distribution functions inside the proton. 
We use the following parameters: ^(lGeV 2 ) = 0.67, ^(lGeV 2 ) = 0.35, ^(lGeV 2 ) = 
0.091, <' p (lGeV 2 ) = 0.029; AT = A d 2 ' p and A d 2 ' n = A^' p , which follow from u n (x,p 2 ) = 
d p (x,p 2 ) and u p (x,p 2 ) = d n (x,p 2 ), respectively |4*H] . 
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Another needed matrix element is [17] 

(p\uu - dd\p) =2M N — = 1.3 GeV . (113) 

The isospin symmetry breaking parameter for the quark condensate is 7 = —0.008 (cf. |12j). 
and the mass parameters of higher resonances are m p i = 1.465 GeV and = 1.649 GeV 
[4~5] . respectively. For the coupling constants we take the values g p7T7T = 6.0, g pl = 5.2 and 
gun = 3 9pi HH- 
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